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I. Introduction 


The general objective of this project is to advance our theoretical understanding of lo's 
atmosphere by studying how various atomic and molecular species are lost from this atmosphere 
and are distributed near the satellite and in the circumplanetary environment of Jupiter. The 
project is divided into well-defined studies described for the likely dominant atmospheric gases 
involving species of the SO 2 family (SO 2 , SO, O 2 , O, S) and for the trace atmospheric gas 
atomic sodium. The relative abundance of the members of the SO 2 family and Na (and its parent 
NaX) at the satellite exobase and their relative spatial densities beyond in the extended corona of 
Io are not well known but will depend upon a number of factors including the upward transport 
rate of gases from below, the velocity distribution and corresponding escape rate of gases at the 
exobase, and the operative magnetospheric/solar-photon driven chemistry for the different gases. 
This question of relative abundance will be studied in this project. 

In order to address this question, we will undertake theoretical modeling studies for the 
distribution and time variability of these exospheric gases in lo's corona/extended clouds and 
thereby evaluate the importance of various physical processes that shape their relative 
abundances and escape rates. Our primary objective will be to study near-Io emission 
observations for O, S, and Na, most of which have already been acquired and some of which are 
scheduled to be acquired in 1996-1997 as part of the larger coordinated International Jupiter 
Watch Observational Campaign to support the Galileo mission. A secondary objective will be to 
continue the study of various larger-spatial-scale groundbased sodium and spacecraft (Voyager 
and Galileo) S02 + observations in order to address related issues and to lay the groundwork for 
larger-spatial-scale O and S observations likely to be obtained in the near future. The proposed 
studies are of scientific importance in understanding (1) the atmosphere of the satellite, (2) the 
interactions of the magnetospheric plasma and the atmosphere, (3) the nature and composition of 
the heavy ion sources for the plasma torus, (4) the impact of these gases on the larger 
magnetosphere, and (5) the spatial distribution of these gases in the magnetosphere and beyond 
in the larger solar wind environment. 

Near-Io observations for this project will be made available in four collaborative efforts 
established with 

(1) F. Scherb of the University of Wisconsin -Madison who from groundbased facilities in 
1990-1996 very successfully obtained synoptic observations of [O I] 6300 A emission 
near Io and in 1996-1998 will continue these synoptic observations, search for [O I] 5577 
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A emission, and add a new Fabry-Perot program element using the new WIYN telescope 
at Kitt Peak to obtain very-high spatial resolution images near lo in the [O I] 6300 A and 
Na 5890 A emission lines, 

(2) G. E. B allester of the University of Michigan who with HST has acquired cycle IV 
and will be acquiring cycle V observations for O and S near Io in various UV emission 
lines, 

(3) L. M. Trafton of the University of Texas who has obtained in 1984-1989, and in his 
ongoing program in 1995 and 1996, groundbased observations for the north-south spatial 
distribution and spectral line shape of sodium (5890 A, 5896 A) emissions near Io, and 

(4) N. M. Schneider of the University of Colorado-Boulder who obtained in 1987 from 
groundbased facilities an extensive set of observations for the east-west spatial 
distribution and spectral line shape of sodium emissions near Io, which exactly overlap 
the October 1987 observations of Trafton. 

These near-Io emissions exhibit time variability with Io System III longitude and Io east-west 
location. The general three year plan of research for these studies is outlined in Table 1. 

II. Summary of Work Performed in the Fourth Bi-Month Periods 

Work accomplished in fourth bi-monthly period includes (1) the completion of an initial 
study for the SO 2 source rate at Io exobase based on presently available Galileo spacecraft data 
and (2) the completion of a study for the total neutral energy loss rate of atomic oxygen and 
sulfur from Io and the implications for local magnetospheric interactions. This research is 
summarized briefly below and is the subject of two papers that are nearing completion and are to 
be submitted for publication in the fifth bi-monthly period. In addition, earlier research work 
included in a paper entitled" Io's Sodium Corona and Spatially Extended Cloud: A Consistent 
Flux Speed Distribution" (Smyth and Combi 1997) appeared in the March 1997 issue of Icarus . 
A copy of this paper is attached to this report. 

To estimate the source strength of SO 2 at the Io's exobase, we have used information for 
the spatial distribution of SC> 2 + obtained from the Galileo spacecraft in its December 7, 1995 
encounter with Io when it flew downstream of the satellite through Io's magnetospheric wake 
with a closest approach distance above the surface of -900 km. An S02 + density profile was 
deduced from magnetic field fluctuations with periods of ~2 to 3 sec measured by the Galileo 
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magnetometer and interpreted as ion cyclotron waves produced by fresh S 02 + Iogenic pickup 
ions created near Io. Calculations for the SC> 2 + density profile along the Galileo trajectory were 
undertaken using the AER neutral cloud model for SO 2 , where S02 + is produced from SO 2 by 
electron impact and charge exchange reactions in the plasma torus. By adopting for SO 2 the 
isotropic incomplete collisional cascade flux source velocity distribution recently deduced for 
atomic sodium at the exobase (Smyth and Combi 1997), we were able to match the minimum 
S0 2 + density profile determined by Huddleston et al. (1997) in their analysis of the ion cyclotron 
waves and the S02 + model calculated density profile for an SO 2 source rate at Io's exobase of -4 
x 10 27 molecules sec"! (225 kg sec" 1) which created in the magnetosphere a total SC> 2 + 
production rate of 1.1 x 10 26 ions sec"! (29 kg sec"!). Most (>90%) of the SO 2 that undergoes 
interactions in the plasma torus is, however, rapidly dissociated primarily by electron impact 
producing O, S, SO, and O 2 which subsequently undergo ionization and charge exchange 
reactions in the plasma torus producing additional mass and energy pickup plasma loading rates 
that are larger than the S02 + production rate of 29 kg sec"! by more than an order of magnitude. 
This substantial exobase source rate for SO 2 of ~4 x 10 27 molecules sec"! j s a factor of only 
about two less than the estimated SO 2 source rate of 8.9 x 10 27 s _1 suggested from the earlier 
analysis of the 0( J D) 6300 A emission brightness near Io by Scherb and Smyth (1993). About 
half of the molecules in the incomplete collisional cascade source distribution for SO 2 are not 
energetic enough to escape the Lagrange sphere of the satellite (-10,550 km or 5.81 Ri 0 ) and 
hence will populate ballistic orbits above the exobase and create a density gradient in the satellite 
corona. This density gradient plays a central role in matching the steepness of the deduced SC> 2 + 
density profile at closer distances to Io (within a radius of -3 satellite radii). The agreement of 
the model and observationally derived SC> 2 + density profiles thus implies that sodium and SO 2 
have a very similar incomplete collisional cascade flux velocity distribution at the satellite 
exobase as would be expected if sodium atoms (a trace species in the atmosphere) were produced 
in the same volume by the same incomplete collisional cascade processes. This agreement also 
implies that the collisional cascade flux velocity distribution at the exobase for other species, 
such as atomic oxygen and atomic sulfur, would also be similar. 

Estimates for the rate at which energy is lost by O and S atoms escaping from near Io's 
exobase because of the interaction of Io's atmosphere with the corotating magnetospheric plasma 
have been made for three source processes: (1) incomplete collisional cascade (2) slow- velocity 
charge exchange and direct ejection, and (3) fast- velocity charge exchange. For the incomplete 
collisional cascade, the energy loss rate is estimated to be 7.5 x 10 9 watts. If this energy loss rate 
is as expected -20% of the total incomplete collisional cascade energy deposition rate in the 
heating layer, then it implies that 56% of the upstream ion kinetic energy flow rate of 6.7 x 10 10 
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watts for a reference cross sectional area of an lo disk area is required for this process. This is a 
factor of 2.8 times larger than the previously adopted value of 20% and implies that the effective 
deflection of magnetospheric plasma out of the interaction region near Io is less than previously 
thought for the collisional cascade process. For slow-velocity charge exchange and direct 
ejection (centered about 20 km/sec), the energy loss rate is estimated to be 3.3 x IO 10 watts, 
which is 50% of the reference upstream ion kinetic energy flow rate. For fast-velocity charge 
exchange (centered about ~60 km/sec), the energy loss rate is much larger and is the dominant 
energy loss mechanism with an estimated value of 6.04 x 10 1 1 watts, which is -900% of the 
reference upstream ion kinetic energy flow rate. The total estimated neutral energy loss rate for 
all three processes (including radiative loss processes) is 6.75 x 10 1 1 watts and is too large by an 
order of magnitude to be provided by the reference upstream initial ion kinetic energy flow rate. 
The required neutral energy loss rate (i.e., power) is, however, matched by the magnetic field 
energy density flow rate of 7.9 x 10 11 watts for the upstream reference area and implies that the 
charge exchange energy loss rates are capable of producing a large pickup current as well as a 
significant reduction in the local planetary magnetic field in the interaction region near Io. The 
magnitude of the current and the reduction in the local magnetic field will depend critically upon 
the volume of the interaction region established by the solution of the three-dimensional 
magnetospheric flow problem past Io including these complex plasma-neutral interactions. 
Rough estimates suggest a pickup current in the range of -4 x 10 6 to 2 x 10 7 A and a reduction 
(AB) in the local magnetic planetary field of -450 nT. This estimated reduction of the magnetic 
field is similar to the remaining and unexplained AB of -400-500 nT determined in a recent 
analysis (Khurana et al. 1997) of the magnetic field depression measured near Io by the Galileo 
Magnetometer (Kivelson et al. 1996 a,b) and attributed by them to an internal magnetic dipole 
field for Io. Hence, we conclude that the reduction of the planetary magnetic field measured by 
Galileo near Io may be a direct reflection of the local charge exchange source and need not 
require an internal magnetic field for the satellite. 
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Table 1 


Three Year Plan of Research for Studies for the Loss of Atomic and Molecular Species from Io 


Subject 


Year 1 


Year 2 


Year 3 


Studies for the SO 2 Family Analyze HST (cycle IV) UV 

data for 0 and S and available 
[O I] 6300 A synoptic data for 
O using the O, S, SO and SO 2 
cloud models; improve model 
execution time, update chemistry 
and refine the model description 
of the plasma torus. 


Analyze the HST (cycle V) UV 
data for O and S; initiate analysis 
of Fabry-Perot image data for 
[01] 6300 A and [01] 5577 A 
(if relevant); re-analyze Voyager 
SOj data; determine sources rates 
and constraints on O, S, SO and 
SO 2 for the individual studies. 


Complete analysis of UV and optical 
data; undertake the comparative and 
collective assessment of the individual 
studies for O, S, Na and SOj to probe 
the nature of the atomic and molecular 
species in Io's atmosphere and their 
implications for the Jupiter system. 


Studies for Sodium 


Continue the analysis of the 1987 
east-west emission data set; initiate 
analysis of the same-date 1987 
north-south emission data; refine 
model description of the plasma 
torus; re-evaluate the Na source 
at Io for the directional feature. 


Complete analysis of the 1987 east- west emission data set; undertake analysis' 
of other select years of the north-south emission data set; determine the nature 
and variability of the Na source conditions and their dependence on east-west 
and System III effects; analyze Fabry-Perot images for sodium and compare 
with [O I] 6300 A images; assess the importance of the electron impact 
excitation and/or nonuniform gas distributions as a cause for asymmetric 
brightness distributions about Io. 
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A data set composed of different groundbased observations 
for lo’s sodium corona and spatially extended sodium cloud 
and covering the spatial range from lo’s nominal exobase of 

I. 4 satellite radii to east- west distances from Io of ± 100 satellite 
radii (jR I o ) is used to investigate the velocity distribution of 
sodium at the exobase. The data set is composed of the novel 
1985 eclipse measurements of Schneider et al. (1991, Astrophys. 

J. 368, 298-315) acquired from —1.4 to —10 R l0 , the 1985 
east-west emission data of Schneider et al. acquired from —4 
to —40 R i0 , and sodium cloud image data acquired near lo’s 
orbital plane from —10 to —100 R lo by a number of different 
observers in the 1976 to 1983 time frame. A one-dimensional 
east-west profile that contains Io is constructed from the data 
set and is analyzed using the sodium cloud model of Smyth 
and Combi (1988, Astrophys. J. Supp. 66, 397-411; 1988, 
Astrophys. J. 328, 888-918). When the directional feature in 
the trailing cloud is either north or south of this east-west line 
(i.e., not at the null condition), an isotropic modified [incom- 
plete ( a = 7/3) coliisional cascade] sputtering flux speed distri- 
bution at the satellite exobase with a peak at 0.5 km sec" 1 
provides an excellent fit to the data set for a sodium source of 
1.7 x 10 26 atoms sec -1 . In particular, the model calculation 
reproduces (1) the essentially symmetric column density distri- 
butions exhibited by the eclipse measurements about Io within 
the Lagrange sphere radius (5.85 R l0 , i.e., the gravitational grasp 
of the satellite), (2) the change in the slope of the column 
density observed just beyond the Lagrange sphere radius in the 
east- west profile of the forward cloud, but not in the trailing 
cloud, and (3) the distinctly different east- west brightness pro- 
files exhibited by the forward and trailing clouds in the emission 
data at the more distant (— ±20-100 R l0 ) portions of the cloud. 
In contrast, the speed dispersion at the exobase for either an 
isotropic Maxwell-Boltzmann flux speed distribution or an iso- 
tropic classical (a = 3) sputtering flux speed distribution (which 
has a higher velocity-tail population than the Maxwell- 
Boltzmann, but not as high as the incomplete coliisional cascade 


sputtering distribution) is shown to be inadequate to fit the 
data set. To fit the enhanced trailing east- west profile observed 
when the directional feature is at the null condition, an addi- 
tional enhanced high-speed (—15-20 km see" 1 ) sodium popula- 
tion is required which is nonisotropically ejected from the satel- 
lite exobase so as to preferentially populate the trailing cloud. 
The need for such a nonisotropic high-speed population of 
sodium has also been recognized in the earlier modeling analysis 
of the directional features (Pilcher et al . , 1984, Astrophys. J. 
287, 427-444), in the more recent lower-velocity component 
required in modeling the sodium zenocorona (Smyth and 
Combi, 1991, J. Geophys. Res. 96, 22711-22727; Flynn et al, 
1992, Icarus 99, 115-130), and in the very recent modeling of 
the directional feature reported by Wilson and Schneider (1995, 
Bull. Am. Astron. Soc. 27, 1154). A complete sodium source 
rate speed distribution function at lo’s exobase from 0-100 km 
sec -1 is then constructed by combining the isotropic modified 
[incomplete (c* = 7/3) coliisional cascade] sputtering flux speed 
distribution, the nonisotropic directional feature (lower-velocity 
zenocorona) source (—15-20 km sec" 1 ) , and the higher-speed 
(—20-100 km sec” 1 ) charge-exchange source required to simu- 
late the sodium zenocorona far from Jupiter. © 1997 Academic Press 


1. INTRODUCTION 

Atomic sodium in the Jupiter system originating from 
a satellite source at Io has been observed in the D 2 (5889.95 
A) and D\ (5895.92 A) emission lines during the past 25 
yr from groundbased facilities. Using an observing slit, 
the sodium emission which is excited by solar resonance 
scattering was first discovered in 1972 by Brown (1974) 
very near Io, where its intensity is brightest [—many tens 
of kiloRaleighs (kR)] and where the sodium density is 
dominated by low-speed (—2 km sec" 1 or less) ballistic 
atom orbits in the satellite “corona.” By occultation of the 
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bright region near Io, image observations (Murcray 1978; 
Murcray and Goody 1978; Matson et al. 1978) were first 
acquired in 1976 and 1977 for fainter (— few to —0.5 kR) 
sodium more distant from Io but still near its circular orbit 
(radius of 5.9 Jupiter radii) about the planet and revealed 
the presence of a predominant “forward cloud” and a less 
spatially extensive “trailing cloud” that moved with the 
satellite. This sodium has been characterized primarily by 
a source of low-speed (—2.6-4 km sec -1 ) atoms that have 
sufficient energy to just escape from Io with an excess 
velocity of only —1 km sec' 1 (or so) and thereby remain 
gravitationally bound to Jupiter fairly near the satellite 
orbit. Additional observations (Pilcher et al. 1984; Gold- 
berg et al. 1984) of even fainter (—1 to 0.2 kR) sodium in 
the early 1980’s revealed a “directional feature” attached 
to Io in the trailing cloud that oscillated north and south 
about the satellite plane with a phase and period deter- 
mined by the Io System III longitude angle. This sodium 
source was characterized by atoms with speeds —20 km 
sec' 1 ejected nonisotropically from the satellite so as to 
populate the trailing cloud and the circumplanetary space 
at larger radial distances beyond Io’s orbit. From earlier 
slit measurements in 1974 (Trafton and Macy 1978), fainter 
(—30 R) sodium emissions well beyond Io’s orbit had been 
observed at a radial distance of —60 planetary radii, while 
from more recent images (Mendillo et al. 1990), very faint 
(—1 R) sodium emissions were observed at radial distances 
of —400-500 planetary radii. Sodium at these larger radial 
distances is called the “magneto-nebula” or “sodium zeno- 
corona” and is thought to be populated primarily by a 
nonisotropic charge-exchange source of high speed (—15- 
100 km sec -1 ) atoms at Io with velocity skewed in the 
forward direction of corotational plasma motion past the 
satellite, and secondarily by a narrow forward sodium jet 
produced by a spatially distributed molecular ion source 
(Wilson and Schneider 1994). Most of this sodium escapes 
the Jupiter system, forms a sodium pause in the sunward 
direction at —2300 planetary radii because of solar radia- 
tion acceleration, and is eventually lost to the solar wind 
by photoionization (Smyth and Combi 1991). 

The observations of sodium emissions on many different 
spatial scales in the Jupiter system thus indicate that its 
atomic source at Io’s exobase must have a wide dispersion 
of speeds. Modeling of these observations has in the past 
been mostly undertaken separately for only one of these 
spatial regions at a time. Although the higher velocity 
dispersions for the sodium zenocorona may be reasonably 
well understood because of its large spatial structure and 
the lack of any significant sodium lifetime impact of the 
magnetosphere, a consistent source for the slower sodium 
in Io’s corona and in the forward and trailing clouds near 
its orbit has not been established. The recent determination 
of the sodium spatial profile in the Io corona obtained 
from the groundbased eclipse data of Schneider (1988; 


Schneider et at. i987, 1991) coupled with earlier emission 
observations, however, now provides a viable observa- 
tional base from which it is possible to pursue the nature 
of this slower sodium. The investigation of a consistent 
exobase sodium source for Io’s corona and the forward 
and trailing clouds near its orbit is therefore undertaken 
in this paper. A consistent flux speed distribution at the 
exobase is determined, and the corresponding sky-plane 
spatial distribution of sodium near Io is presented. Sodium 
source information obtained from previous modeling anal- 
ysis of Io’s corona and the forward and trailing clouds is 
first summarized in Section 2. The observational data base 
to be investigated in this paper is presented in Section 3. 
Modeling of an east-west spatial profile determined from 
this observational data base is undertaken in Section 4. 
Discussion and conclusions are presented in Section 5. 

2. EARLIER SODIUM MODELING 

The major modeling analysis studies for the spatial distri- 
bution of sodium near Io and its orbit are summarized in 
Table I. The summary is divided into three observed spatial 
regions: (1) the Io corona located within the satellite La- 
grange sphere (average radius of 5.81 R Io or —3 arcsec), 
(2) the neutral cloud located beyond the Lagrange sphere 
and near Io’s orbit, and (3) the north-south oscillating 
directional feature, observed to trail Io in its orbit. Model- 
ing analysis studies for the Io corona are further subdivided 
into early observations of the average intensity in an 8 X 
3 arcsec slit centered on Io reported by Bergstralh et al. 
(1975, 1977) that indicated an east-west intensity asymme- 
try of —1.25 and later observations for one-dimensional 
column density profiles within the Lagrange sphere re- 
ported by Schneider et al. (1987, 1991). 

2.1. Corona: East-West Intensity Asymmetry 

In Table I, the early studies of Smyth (1983) for sodium 
atoms ejected monoenergetically from Io’s exobase estab- 
lished that small scale structures in the D-line intensity 
profile observed as a function of the Io geocentric phase 
angle (Bergstralh etal. 1975, 1977) could arise from modu- 
lation of the atoms’ escape rate from Io caused by the 
action of solar radiation acceleration in the D-lines. These 
modulations occur primarily for exobase speeds near 2.0 
and 2.1 km sec" 1 , which are near the escape-speed thresh- 
old of the Lagrange sphere. Later studies of Smyth and 
Combi (1987a) showed that the main reason for the east- 
west intensity asymmetry was, however, an east-west elec- 
tric field which altered the plasma properties at Io’s orbit 
so as to increase the sodium lifetime and hence sodium 
abundance when Io was preferentially east of Jupiter. More 
complex modeling studies of Smyth and Combi (1988b) 
constrained the flux velocity dispersion for sodium at Io’s 
exobase by simultaneously fitting the average east-west 



TABLE I 

Summary of Modeling Studies for the Spatial Distribution of Sodium Near Io and Its Orbit 


Spatial Region 

Topic Studied 




Sodium Source 




Sodium Lifetime 

Orbital Dynamics 


Mono-energetic 
Atom Ejection 

Speed 
(km s'*) 

Maxwell-Boltzman 
Flux Distribution 

T V peak 

<K) (kms'*) 

a 

Cascade Flux Speed Distribution 

v b v peak ^max 

(kms'*) (kms'*) (kms'*) 

Angular 

Nature 3 

Plasma Torus 
Description 

Lifetime 

(hrs) 

lo’s Mass 
Included 

Jupiter’s 

Mass 

Included 

Radiation 

Pressure 

Included Reference 

l. Corona 

e-w intensity asymmetry* 1 

2. 0-2.6 

- 

- 

- 

- 


- 

isotropic 

cut-off 

20 

yes 

yes 

yes 

Smyth 1983 


e-w intensity asymmetry* 5 

2.6 

- 

- 

- 

- 



isotropic 

2-D 

variable 

yes 

yes 

yes 

Smyth and Combi 1987a 


e-w intensity asymmetry* 5 

- 

460; 12,300 

0.71; 3.65 

7/3 

<0.4; 2.2 

<0.5; 2.9 

46.6 

band, isotropic 

2-D 

variable 

yes 

yes 

yes 

Smyth and Combi 1988b 


column-density profile 0 

2,6 

- 

- 

- 

- 



isotropic 

2-D 

varialbe 

yes 

yes 

yes 

Smyth and Combi 1987b,c 


column-density profile 0 

- 

- 

- 

3 

0 d 


<50 

isotropic 

cut-off 

3,~ 

yes 

no 

no 

McGrath 1988 


column-density profile 0 

- 

1000 

1.04 

- 

- 



isotropic 

uniform 

oo 

yes 

no 

no 

Summers et al. 1989 


column-density profile 0 

3.0 

- 


- 

- 



isotropic 

uniform 

? 

yes 

no 

no 

Ip 1990 


column-density profile 6 

- 

1500 

1.28 

- 

- 


- 

isotropic 

uniform 

oo 

yes 

no 

no 

Schneider et al. 1991 

II. Cloud 

general spatial nature 

3.5 

- 

. 

- 

- 


- 

isotropic 

uniform 

30, 47 

yes 

yes 

no 

Carlson et al. 1975 


general spatial nature 

- 

5000 

2.3 

- 

- 



isotropic 

long-lived limit 

- 

no 

yes 

no 

Fang et al. 1976 


evolution and 2-D nature 

2,3 

- 


- 



- 

isotropic 

cut-off 

50 

yes 

yes 

no 

Smyth and McElroy 1977 


I-D brightness morphology*^ 

- 

- 

- 

3 

4 

4 

oo 

I-L hemispheres 

uniform 

28 

yes 

yes 

no 

Matson et al. 1978 


2-D brightness morphology* 1 

2.6,3 

- 

- 


- 


- 

I, I-T hemisphere* 

cm-off 

15,20 

yes 

yes 

no 

Smyth and McElroy 1978 


e-w orbital asymmetry* 

2.6 

- 

- 

- 

- 



I-T hemisphere 1 

cut-off 

20 

yes 

yes 

yes 

Smyth 1979 


radial and vertical structure 

(3,5, 7,9,1 1,1 3) k 

- 

- 

- 

- 


- 

I hemisphere 

uniform 

56 

no 

yes 

no 

Macy and Trafton 1980 


1-D brightness morphology 

- 

- 

- 

3* 

4* 

4* 

<50 i 

I-L hemisphere* 

1-D 

variable 

yes 

yes 

no 

Goldberg et al. 1980 


e-w orbital asymmetry 

2.6 

- 

- 

- 

- 

- 


1 hemisphere 

cut-off 

20 

yes 

yes 

yes 

Smyth 1983 


2-D brightness morphology 

2.6 

- 

- 

- 

- 

- 


band, isotropic 

2-D 

variable 

yes 

yes 

yes 

Smyth and Combi 1988b 

HI. Directional Feature 

spacetime structure 01 

~20 


. 


. 

. 

- 

O, -± to lo’s motion 

2-D 

variable 

yes 

yes 

yes 

Pilcher etal. 1984 


collision cross sections 

- 

- 

- 

3 

0° 



“JL to lo's motion 

- 

- 

- 

- 

- 

Sieveka and Johnson 1984 


- nearer zenocorona structure 

- 

- 

~2G± 12P 

- 

- 



tangential ± isotropic 

photoionization 

-400 hr 

no 

no; Sun yes 

yes 

Smyth and Combi 1991 


nearer zenocorona structure 

- 

- 

“20± 12P 

- 

- 

- 

- 

tangential ± isotropic 

photoionization 

-400 hr 

no 

yes 

no 

Flynn et al. 1992 


spacetime structure 

- 

- 

20± 10-209 

- 

- 

- 


tangential ± isotropic 

none 

oo 

no 

yes 

no 

Wilson and Schneider 1995 


a regarding Jo's exobase, L=leading, I=inner, T=trailing, O=outer. 
b east- west intensity asymmetry data of Bergstralh et ai. {1975, 1977). 
c early (i.e., partial) Na eclipse data set from Schneider et al. (1987). 

d used cut-off energy: Emin < E < «; Emin = k T eX obase! T exo base = 1500 K (i.e. a 1.04 km s'* cut-off speed). 
e complete eclipse data set from Schneider (1988) and Schneider et al. (1991). 

^ two sodium cloud images of Matson et al. (1978). 

8 I-L hemisphere centered 30° longitude (0° longitude facing Jupiter, 90° longitude is the leading point in the orbit). 
** fifty-six sodium cloud images of Murcray (1978) and Murcray and Goody (1978). 


1 I-T hemisphere centered on -40° longitude. 

■* east-west orbital asymmetry data of Goldberg et al. ( 1978). 
k velocity components equally weighted. 

* distribution parameters from Carlson (1995, private communication); I-L hemisphere centered on 45° longitude. 
m from images of Pilcher et al. (1984). 

n used cut-off energy of 0.5 ev (i.e., Na cut-off speed of 2.0 km s'* ). 

P tangential speed to lo's orbit at Io position ± most probable speed of an isotropic Maxweli-Boltzmann. 

9 tangential speed to lo’s orbit at fo position ± most probable speed of an isotropic Gaussian. 
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intensity asymmetry and also the general spatial morphol- 
ogy of the forward sodium cloud, located on a much larger 
spatial scale well beyond the Lagrange sphere. These stud- 
ies showed that the sodium ejection speed at the exobase 
required to fit the east-west intensity asymmetry is double- 
valued, having a lower value of ^1 km sec -1 and a higher 
value in the range 2.6-3.65 km sec -1 . For a Maxwell- 
Boltzmann flux distribution, the lower and higher most 
probable speed values were 0.71 km sec -1 (T = 460 K) 
and 3.65 km sec -1 (T = 12,300 K). Neither distribution 
was, however, suitable for properly populating the forward 
cloud. The lower value produces essentially only ballistic 
atom orbits which could not populate the forward cloud, 
while the higher value was significantly larger than the 
nominal ~2.6 km sec -1 characteristic monoenergetic veloc- 
ity required to reproduce the proper spatial morphology 
of the forward cloud as a function of the Io geocentric 
phase angle. For a Maxwell-Boltzmann flux distribution 
with a more nominal thermal exobase temperature in the 
range —1000-2000 K, the calculated east-west intensity 
ratio was much higher than the observed value, with the 
atoms still contributing primarily to the corona density 
and again far too deficient in energy to contribute any 
significant sodium to the forward cloud. For the preferred 
(a = 7/3) modified-sputtering distribution of Smyth and 
Combi (1988b) with source strength —2 X 10 26 atoms sec -1 , 
the lower and higher speed values were <0.5 and —2.9 km 
sec -1 , respectively, with the latter value being preferred 
because of its closer proximity to the —2.6 km sec -1 charac- 
teristic velocity for the forward cloud. Interestingly, how- 
ever, it is actually the lower value that will be shown in 
this paper to reproduce the correct spatial profile for so- 
dium both within the Lagrange sphere and beyond in the 
more distant neutral cloud. 

2.2. Corona: Column Density Profile 

In Table I, modeling studies of Smyth and Combi 
(1987b, c) determined that typical forward cloud brightness 
data for the sodium cloud could be properly simulated 
well beyond the Lagrange sphere radius of —5.81 R io by 
a sodium source of X 10 26 atoms sec -1 ejected monoen- 
ergetically from Io’s exobase with a characteristic velocity 
of ~2.6 km sec -1 . They also established that this same 
sodium source reproduced the column density profile of 
Schneider et al (1987) within the Lagrange sphere down 
to a radius of —3.5 R\ 0 . For a radius smaller than '-3.5 R Xo , 
the calculated profile was lower than the observed profile, 
indicating that lower (ballistic) velocity components are 
required, in addition, as part of a more realistic flux velocity 
dispersion. A similar behavior for the simulated column 
density profile, with an even more dramatic departure from 
the observed profile both inside and outside the Lagrange 
sphere, was also later shown by a model calculation of Ip 


(1990), who assumed an exobase speed of 3 km sec -1 but 
did not include the gravity of Jupiter so as to properly 
include the near zero escape speed conditions for sodium 
at the Lagrange sphere. Adopting for sodium atoms at the 
exobase a simple (i.e., binding-energy velocity = 0) 
classical sputtering energy distribution with a low energy 
cut-off and also excluding Jupiter’s gravity, McGrath 
(1988) modeled the column density within the Lagrange 
sphere and produced a profile with a slope slightly less 
steep than the observation for an infinite sodium lifetime 
and a slope somewhat steeper than the observation for a 
sodium lifetime of 3 hr. Alternatively adopting a Maxwell- 
Boltzmann flux distribution, assuming an infinite sodium 
lifetime, and similarly excluding Jupiter’s gravity. Summers 
etal. (1989) and Schneider etal. (1991) modeled the column 
density within the Lagrange sphere region and produced a 
profile that reasonably well matched the observed profile 
for an exobase temperature, respectively, of 1000 K based 
on the partial eclipse data set (Schneider et ai 1987) and of 
1500 K based on the complete eclipse data set (Schneider et 
al 1991). Although these different flux velocity distributions 
reasonably fit the observations within the Lagrange sphere, 
it is clear from the earlier studies of Smyth and Combi 
(1988b) that these Maxwell-Boltzmann distributions are 
energetically deficient and inappropriate for populating the 
neutral cloud and, furthermore, that the more energetically 
promising sputtering distribution cannot be investigated ad- 
equately near or beyond the Lagrange sphere radius without 
properly including the gravity of Jupiter, solar radiation ac- 
celeration, and the spacetime variable sodium lifetime in the 
plasma torus. This study will be undertaken in Section 4. 

2.3. Sodium Cloud 

The early studies in Table I for the sodium cloud were 
general in nature, probing its poorly documented spatial 
and angular extent about the planet. Based upon the solar 
resonance scattering excitation mechanism for sodium 
(Bergstralh et ai 1975) and limited angular extent data 
determined by slit-averaged intensity data, Carlson et al. 
(1975) undertook monoenergetic (3.5 km sec -1 ) model cal- 
culations and estimated that the sodium cloud lifetime 
(assumed to be spatially uniform) was likely determined 
by electron impact ionization by the (then very poorly 
characterized) plasma in the planetary magnetosphere. 
This general picture for the cloud was confirmed by more 
extensive model calculations performed by Fang et ai 
(1976) and Smyth and McElroy (1977), the latter of which 
explored the time evolution and two-dimensional nature 
of the cloud for exobase velocities near the Io Lagrange 
escape speed. The acquisition of sodium cloud images in 
late 1976 and early 1977 brought this subject into dramatic 
focus. For a classical sputtering flux distribution that 
peaked at 4 km sec -1 , Matson et al. (1978) successfully 
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modeled a one-dimensional east-west brightness profile 
(derived from a cloud image) which extended from Io in 
the forward cloud to —80 R lo and in the trailing cloud to 
—40 R\ Q but which excluded sodium emission within Io’s 
corona. The analysis (Smyth and McElroy 1978) of a much 
larger sodium cloud image data set (Murcray 1978) also 
indicated that the forward cloud could be characterized 
by an exobase ejection speed of —2.6 km sec” 1 and that 
its changing intensity pattern could be understood as the 
changing viewing perspective of an approximately steady 
state cloud on the sky plane as Io moved on its orbit around 
Jupiter. The observed predominance of the forward so- 
dium cloud over the trailing cloud was accomplished in all 
these models by limiting the exobase source area to a 
hemisphere (see Table I) and by limiting the assumed 
spatially uniform lifetime so as to dynamically select so- 
dium atom orbits that would primarily populate the for- 
ward cloud. Additional modeling studies by Macy and 
Trafton (1980) of the radial and vertical cloud structure 
on a larger spatial scale indicated source dispersion speeds 
at least up to 13 km sec” 1 were, however, required to 
explain a variety of other observations. Additional model 
studies (Smyth 1979, 1983) showed that the newly discov- 
ered east-west orbital asymmetry of the sodium cloud 
(Goldberg et al 1978) was not source related but was due 
to the perturbing action of solar radiation acceleration 
on the sodium atom orbits. Adopting a one-dimensional 
radially dependent sodium lifetime in the plasma torus 
based upon limited Voyager spacecraft data and an asym- 
metric exobase source for a classical sputtering distribution 
with a peak velocity at 4 km sec" 1 , Goldberg et al. (1980) 
successfully modeled a one-dimensional east-west bright- 
ness profile acquired during the Voyager 1 encounter for 
distances extending from Io in the forward cloud to —80 
Rio and in the trailing cloud to —30 R Xo . Later modeling 
by Smyth and Combi (1988b) using a more accurate two- 
dimensional space- and time-dependent sodium lifetime 
in the plasma torus and an isotropic (or near isotropic) 
exobase sodium source of —2 X 10 26 atoms sec" 1 demon- 
strated that the predominant forward cloud was caused by 
the highly radially dependent sink for sodium in the plasma 
torus and not by a nonisotropic source. The deduced char- 
acteristic or most probable exobase speed for the more 
definitive modeling of the forward sodium cloud above is, 
therefore, in the range —2.6 to 4 km sec” 1 and is much 
larger than required to characterize the sodium column 
density profile in Io’s corona. A new flux speed distribution 
is therefore needed for consistency and is determined in 
Section 4. 

2.4 . Directional Feature 

In Table I, observations acquired in 1980 and 1981 by 
Pilcher et al. (1984) for weaker D-line emissions in the 


trailing portion of the sodium cloud allowed them to dis- 
cover an elongated feature in the brightness distribution 
that on the sky plane was directed away from Jupiter and 
was inclined sometimes to the north and sometimes to 
the south of the satellite’s orbital plane. The north-south 
direction of the feature was shown to be correlated with 
Io’s magnetic longitude and suggested a formation mecha- 
nism involving the oscillating plasma torus. Modeling anal- 
ysis by Pilcher et al. indicated that the feature resulted 
from a high-velocity (—20 km sec” 1 ) sodium source that 
was at near right angles to Io’s orbital motion with a source 
strength required on the outer satellite hemisphere of 
— 1 X 10 26 atoms sec” 1 . This peculiar directionality of the 
source was investigated by Sieveka and Johnson (1984), 
who concluded that it was likely produced by direct colli- 
sional ejection of neutral sodium from the exosphere by 
the corotating plasma flow past Io. Modeling of the sodium 
zenocorona (Smyth and Combi 1991; Flynn et al 1992) 
showed that it was consistent with a two-component exo- 
base source: a similar high-velocity (—20 km sec" 1 ) sodium 
source of —1 X 10 26 atoms sec" 1 for the spatial distribution 
nearer the planet and an even higher-velocity (—57 km 
sec" 1 ) sodium source of —2 X 10 26 atoms sec" 1 for the 
spatial distribution further from the planet. Both source 
components, however, were based on ion-neutral charge 
exchange processes in Io’s exosphere and were hence com- 
posed of a speed tangential to Io’s orbit at Io’s position 
plus an isotropic Maxwell-Boltzmann distribution with a 
most probable speed of about one-third of the tangential 
speed in the Jupiter frame (i.e., one-third of —37 and —74 
km sec" 1 , respectively). In Table I, the lower component 
is therefore symbolically denoted in the Io frame by 20 ±12 
km sec" 1 . Recent modeling of the directional features has 
also been reported by Wilson and Schneider (1995), who 
used a similar lower component source denoted in Table 
I by 20 ± 10-20 km sec" 1 , where the isotropic portion of 
their source may be variable in magnitude. A spatially 
distributed molecular-ion source has also been used to 
model the sodium zenocorona (Flynn 1993), but appears 
to be of secondary importance due to its smaller source 
rate and erratic presence. 

3. OBSERVATIONAL DATA BASE FOR MODELING 

To describe the entire spatial distribution of sodium in 
Io’s corona and beyond in the extended neutral clouds, 
three different types of sodium observations obtained on 
very different spatial scales are combined. For Iocentric 
distances, the combined data set is composed of the novel 
1985 eclipse measurements of Schneider et al. (1991) ac- 
quired from —1.4 to —10 R Xo , the 1985 east-west emission 
data of Schneider et al (1991) acquired from —4 to —40 
Rio , and sodium cloud image data acquired near Io’s orbital 
plane from —10 to —100 R Xo by a number of different 
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TABLE II 

Io Eclipse and Emission Observations for 1985 


Date 

UT 

Midpoint 

Io Geocentric 
Phase Angle 
Range 
(deg) 

Io System III 
Longitude 
Range 
(deg) 

Spectrum ID 

Type of 
Observation 

Eclipse Emission 

Dominant 

Spatial 

Profile 

Directional 

Feature 

Orientation 

Enhance 

Doppler 

Signature 

D 2 Emission Profile Power Law Fit ^ 
Exponent Amplitude (kR) 

East West East West 

Profile Profile Profile Profile 

August 27 

0714 

61.4 ± 0.7 

29.9 ± 2.3 

85gl88 

1 

forward 

null 

trailing/forward 

1.67(F) 

1.57 (T) 

191 

124 


0720 

62.3 ± 11.0 

32.9 ±35.8 

a 

1 

— 

— 

— 






0830 

72.2 ± 1.1 

65.1 ± 3.4 

85gl96 

2 

symmetric 

south 

no 

1.85 (S) 

1.80 (S) 

169 

142 

September 13 

0641 

276.6 ± 0.7 

194.7+ 2.3 

85h032 

3 

trailing 

null/north 

trailing 

... 

- 

... 

... 

September 14 

0245 

87.7 ± 0.7 

31.4 ± 2.3 

85 h 102 

4 

trailing 

null 

trailing 

1.23 (S/T) 

1.80 (S/F) 

89 

188 


0326 

93.5 ± 4.2 

50.5 ± 13.7 

b 

2 

— 

— 

— 






0416 

100.6 ± 0.4 

73.6 ± 1.4 

85hl 13 

5 

g 

south 

no 

1.27 (T) 

... 

66 

... 

September 15 

0316 

294.7 ± 0.7 

353.7 ± 3.0 

85hl52 

6 

trailing 

null 

trailing 



... 

... 

... 


0500 

309.4 ± 11.5 

41.7 ±37.4 

c 

3 


... 

"* 





September 21 

0604 

100.5+ 1.6 

112.5+ 5.2 

d 

4 

- 

... 

... 





September 23 

0230 

11 7.2 ± 0.7 

267.6 ± 2.3 

85h433 

7 

forward 

north 

no 

2.05 (T) 

1.57 (F> 

342 

135 


0301 

121.6 ± 0.4 

281.9+ 1.2 

85h436 

8 

forward 

north 

no 

1.96 (T) 

1.54 (F) 

283 

138 


0356 

129.3 ± 2.0 

307.3 ± 6.7 

e 

5 

— 

— 

— 






0534 

143.1 ± 0.7 

352.9 ± 2.3 

85h457 

9 

forward 

null 

no 

2.16 (T) 

1.64(F) 

374 

165 


a. Eclipse 1 : 85gl79, 85g181, 85gl85, 85g188. 85g192, 85g193, 85gl96 

b. Eclipse 2 : 85h103, 85h104, 85h105, 85h106, 85M 07, 85htOQ, 85h109 r 85h1 10. 85hl 12 

c. Eclipse 3 : 85h153, 85h154. 85H155, 85h157, 85h159, 85h162, 85h163 

d. Eclipse 4 : 856287, 856288, 85h289, 85h290. 85h291, 856292, 856293. 85h294, 85h295 

e. Eclipse 5 : 856441, 856442, 856443, 85h444, 85h445, 856446. 856447, 85h448. 856449, 856450 

f. Profile points inside of 4 Rj 0 are excluded; power law fit A r'P, where A is the amplitude, (3 is the exponent, and r is in units of a* F - forward cloud; S = symmetric turning point; T = trailing cloud. 

g. Not sufficient data west of Io to compare spatial profils (see Table HI). 

(Note one lower-bound data point from eclipse 4 at a distance from the center oflo of 1.17 Ri q is excluded in the analysis since it is well within the nominal exohase radius of 1 .4. R] 0 ) 


observers in the 1976 to 1983 time frame. For modeling 
purposes in Section 4, a one-dimensional east -west profile 
centered on Io is constructed from this data set. From the 
data of Schneider et al. (1991), five higher quality eclipse 
profiles and nine higher quality emission profiles have been 
selected, and their observational dates, times, Io angular 
parameters, spectral ID numbers, and the numbering of 
these interleaved observations as adopted in this paper 
are summarized in Table II. The D 2 brightnesses for the 
emission profiles in Table II, previously published only in 
a graphical format, are given numerically in Table III as 
provided by Schneider (1990, 1995, both private communi- 
cations). For the sodium cloud data, fourteen images in the 
D 2 emission line acquired in the 1976-1983 time interval 
(Murcray 1978; Murcray and Goody 1978; Matson et al 
1978; Goldberg et al 1980, 1984; Morgan 1984, private 
communication) have been selected, with values for the 
east-west D 2 brightness profiles of the forward and trailing 
clouds extracted and summarized in Table IV. 

The eclipse observations provide the most accurate in- 
formation in the radial interval from —1.4 to 6 R io for the 
atomic sodium column density profile in the corona within 
the Lagrange sphere of Io ( i.e., a radius of 5.81 R {G ) and 
yield an essentially symmetric column density profile about 
Io with a power law fit N (1.4 < r < 5.85) = 2.55 X 
10 i2 r-2- 48 7 where N is in units of atoms cm -2 and r is the 
distance from the Io’s center in units of R Xo . This power 
law fit, however, undercuts the eclipse data beyond the 
Lagrange radius, and this reduced slope will later be seen 
to be caused by the dominant planetary gravitational field 


beyond Io’s Lagrange radius. The emission observations 
provide accurate information for the sodium D 2 brightness 
from just within the Lagrange sphere outward into the 
nearer portion of the sodium cloud (i.e., —4-40 i? l0 from 
Io’s center) with a power law fit 7 D2 (r>4) = 101 r“ L45 , 
where l Di is in units of kiloRayleighs (kR) and where the 
D 2 brightness of —100 kR as r approaches Io’s surface is 
consistent (see Brown and Yung 1976) with the maximum 
sodium column density of —1 X 10 12 atoms cm -2 deduced 
from the eclipse data. The brightnesses for the different 
observed profiles in Table III, however, vary by a factor 
of —3 to almost 5 at the same distance from Io but have 
error bars that are no larger than ±30%, suggesting that 
the large variation is real and likely correlated with Io 
geocentric phase angle, Io System III longitude, and the 
east-west asymmetry in the plasma torus, as is the case 
for sodium cloud image data. The closer spatial regions 
covered by the eclipse data and the near Io emission data 
are, however, masked in the cloud images by a circular 
(or nearly circular) occulting mask of —10 R io in radius 
centered on Io. The sodium cloud image data in Table IV, 
as illustrated in Fig. 1, hence contribute to the east- west 
brightness profiles for Iocentric distances from —10 to 
—100 7 ? l0 . 

The structure of the sodium cloud emission brightness 
on the sky plane has been historically divided into a for- 
ward cloud, so-called because it appears ahead of the satel- 
lite in its orbit (i.e., in Fig. 1 located right (west) of Io in 
image A and left (east) of Io in images B and C), and 
a corresponding trailing cloud that appears behind the 
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TABLE III 

Emission Data for 1985 


Radial 
Distance 
From Io 
(satellite radii) 

D 2 Intensity$(kR) 

Emission 1 
61.4° 

27 August 
(85gl88) 

Emission 2 
72.2° 

27 August 
(85gl96) 

Emission 4 
87.7° 

14 September 
(85hl02) 

Emission 5 
100.6° 

14 September 
(85hl 13) 

Emission 7 
117.2° 

23 September 
(85h433) 

Emission 8 
121.6° 

23 September 
(85h436) 

Emission 9 
143.1° 

23 September 
(85h457) 

Emission 3* 
276.6° 

13 September 
(85h032) 

Emission 6 
294.7° 

15 September 
(85hl52) 


-29.92 





0.6910.17 

0.9310.18 

0.6610.16 




-29.65 









1.1210.19 


-24.48 








1.4010.21 



-21.76 





1.1910.19 

0.8410.17 

0.9510.18 




-16.32 





1.3810.20 

2.1010.26 


2.2210.27 



-15.78 







1.7410.23 


2.1510.26 


-14.96 




1.91 ±0.24 







-10.88 

3.16 ±0.35 




2.9010.33 

3.2910.36 

3.3810.37 

4.4010.46 


T 

-9.52 



3.1710.35 







West 

-6.80 

5.21 ±0.54 

4.5010.47 

6.41 ±0.66 


7.4910.76 

7.7610.79 


11.5411.16 



-4.08 

14.65 ± 1.47 

11.291 1.14 

14.71 1 1.48 




7.5410.77 

28.6912.87 



-3.81 









18.7512.63 


-2.18 

24.31 ±3.40 

18.3912.57 

24.1213.38 


25.6113.59 

25.01 ±3.50 

21.0012.94 

45.4716.37 



-1.09 

27.46 ±3.84 

19.9912.80 

27.1413.80 


29.9514.19 

28.7214.02 

23.8913.34 

48.6116.81 



-0.82 









24.8513.48 

Io 

0.00 

27.44 ±3.84 

19.1912.69 

27.1413.80 


27.9013.91 

30.1614.22 

23.9713.36 

45.5616.38 



0.27 









25.27 1 3.54 


1.09 

25.8713.62 

19.3512.71 

26.6513.73 

23.0313.22 

26.31 ±3.68 

26.39 1 3.69 

22.8413.20 

39.45 1 5.52 



1.36 









22.9013.21 


2.18 

25.1913.53 

16.8712.36 

25.6813.60 

19.6012.74 

22.4213.14 

23.91 ±3.35 

20.2112.83 

33.6014.70 



3.26 









11.5311.61 


4.08 

17.6011.77 

11.921 1.20 

18.101 1.82 

11.8711.20 

15.5411.56 

14.831 1.49 

12.5211.26 

20.0612.01 



5.98 









4.7910.50 

East 

6.80 

8.31 ±0.84 

5.5210.57 

8.1710.83 

5.52 ± 0.57 

7.0910.72 

7.1010.73 

6.0210.62 

7.5210.77 


1 

10.06 









1.6610.22 


10.88 

3.4010.37 

1.9410.25 

3.8410.41 

3.2! ±0.35 

2.47 1 0.29 

2.47 ± 0.29 

2.2610.27 

2.7310.31 



15.50 









0.7710.17 


16.32 



2.8410.32 

1.8110.24 

1.0310.18 

1.1110.19 

0.8210.17 

1.3910.20 



21.76 



2.2510.27 


0.6810.16 

0.7510.17 

0.5210.16 




-23.66 









0.5310.16 


24.48 




1.2410.19 




1.0610.18 



27.20 



1.3910.21 








35.36 



1.2710.20 





1.0010.18 



46.24 








0.8310.17 



^Nonuniform spatial coverage occurs because of different distance intervals adopted to obtain good average brightness values (given different signal to noise ratios) and because 
of signal drop-out associated with constraints imposed on positioning the slit profile on the CCD detector during interleaved eclipse and emission measurements. 

*Calibration uncertain 


TABLE IV 

East-West D 2 Brightness Profiles for Sodium Cloud Image Data 


Io Geocentric Io System III 
Image ID Phase Angie Longitude 

UT Date UT Time Number (deg) (deg) 


East-West Distance from Io for Specified D 2 Brightness Level (satellite radii) 


Image Data Set Reference: Murcray (1978) 



0.5 kR 


1.0 kR 


1.5 kR 


2.0 kR 






Forward 

Trailing 

Forward 

Trailing 

Forward 

Trailing 

Forward 

Trailing 

1976 Nov. 16 0806 

ES 328B 

256 

262 

- 

- 

51 

29 

46 

24 

28 

14 

1977 Jan 27 0024 

ES 3 69 A 

86 

193 

60 

38 

50 

22 

44 

19 

32 

16 

0217 

ES 370D 

102 

245 

>83 

30 

56 

25 

46 

21 

24 

19 

Image Data Set Reference: Goidberg(1988f) 



0.2 kR 


0.5 kR 


I.OkR 


2.0 kR 






Forward 

Trailing 

Forward 

Trailing 


Trailing 

Forward 

Trailing 

1981 May 5 0819 

SIP 418/31-33 

102 

300 

78 

69 

63 

40-66 

45 

37 

24 

21 

May 12 0848 

SIP 420/30-32 

91 

302 

66 

86 

41 

37 

35 

29 

20 

21 

May 13 0346 

SIP 421/21-23 

253 

108 

74-103 

73 

70 

41 

41 

27 

29 

20 

0555 

SIP 421/32-33 

27! 

168 

124 

65 

112 

44 

51 

30 

26 

19 

June 6 0436 

SIP 424/10-12 

103 

300 

81 

75 

68 

71 

36 

26 

20 

23 

Image Data Set Reference: Morgan (I984f) 



0.3 kR 


0.6 kR 


0.9 kR 


1.8 kR 






Forward 

Trailing 

Forward 

Trailing 

Forward 

Trailing 

Forward 

Trailing 

1983 June 13 0714 

i 8492 

274 

230 


39 

58 ■ 

26 

37 

19 

20 

- 

0722 

i 8494 

275 

233 

_ 

35 

52 

22 

40 

20 

20 


0729 

i 8496 

276 

237 

>93 

36 

52 

22 

37 

17 

22 

-- 

0827 

i 8501 

284 

264 

>93 

39 

61 

25 

42 

19 

23 

- 

0949 

i 8509 

296 

301 

>93 

47 

63 

26 

44 

17 

24 

- 

1010 

i 851 1 

299 

311 

>93 

41 

63 

23 

44 

16 

19 

- 


tlmage observational data obtained by private communication. 
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10 SODIUM CLOUD 1981 



FIG. 1. lo sodium cloud images. Three calibrated D 2 emission images of the Io sodium cloud from the JPL Table Mountain Data Set are shown 
to proper scale with Jupiter and Io’s orbit as viewed from earth in 1981 (Smyth and Goldberg 1993). The Io System III longitude and corresponding 
orientation of the trailing directional feature in image A are 247° and north, in image B are 104° and south, and in image C are 178° and only very 
slightly north. An east-west spatial scale of ±100 satellite radii about Io is shown for reference, with tick marks also located at ±50 satellite radii. 
Contour levels for the D 2 brightness, from outside to inside, are 0.2, 0.5, 1, 2, 4, 6, 8, and 10 kR. An occulting mask of ~10 Ri 0 in radius is centered 
on Io so that brightness values within this distance are not accurate. 


satellite. The forward cloud changes in length, brightness, 
and east-west orientation relative to the satellite’s location 
as Io moves about Jupiter (i.e., as a function of the Io 
geocentric phase angle). This change in orientation is well 
documented (Murcray 1978; Murcray and Goody 1978; 
Goldberg et al 1984) and is due primarily to the projection 
upon the two-dimensional sky plane of a three-dimensional 
cloud (Smyth and McElroy 1978), slightly altered from 
mirror symmetry by solar radiation pressure (Smyth 1979, 
1983), which passes through an east symmetric turning 
point between about 65° and 85° to 90° and a west symmet- 
ric turning point at about 235° (Goldberg et al 1984). A 
detailed examination of the emission data of Schneider et 
al (1991), summarized in the last seven columns of Table 
II, shows that the forward and trailing profiles for the 
emission data are quite consistent with this known behavior 
of the Io sodium cloud images. In the trailing cloud, the 
time-dependent change in the north-south inclination of 
the fainter directional feature (see Fig. 1) has been shown 
to be correlated with the System III longitude of Io (Pilcher 
et al 1984; Goldberg et al 1984) with the directional feature 
changing from a south to north inclination (a first null 
point) at an Io System III longitude near 165° and changing 
from a north to south inclination (a second null point) for 


a rather poorly defined Io System III longitude somewhere 
between about 320° and 25°. When the directional feature 
is near the null location, as illustrated in Fig. 1, an increase 
in both the spatial extension and brightening of the trailing 
cloud along the east- west oriented (dashed) line is readily 
apparent. In addition, since the trailing cloud is associated 
with a high-speed Io sodium source (—15-20 km sec" 1 ), 
an increase in the Doppler width of the spectral line in 
the trailing cloud brightness along an east-west slit is also 
expected near the null location and is indeed observed in 
the emission data of Schneider et al (1991), as indicated 
in Table II. 

Information for the sodium D 2 brightness profiles in 
both the forward and trailing clouds is presented in Fig. 
2. The forward and trailing cloud orientation depicted in 
Fig. 2 is chosen for Io near eastern elongation in order 
to facilitate the comparison with the 1985 emission data 
profiles, mostly acquired for Io east of Jupiter. The five 
emission profiles of Schneider et al (1991) for Io east of 
Jupiter are shown by different symbol together with their 
power-law fits (/ D2 (r ^ 4) = Ar~ $) given in Table II. For 
the Io sodium cloud images, the extracted east-west D 2 
brightness profiles in Table IV for the forward and trailing 
clouds are shown by shaded areas, which represent appro- 
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TRAILING CLOUD 


FORWARD CLOUD 



Distance East of lo (R Jq ) Distance West of lo (R jo ) 

FIG. 2. East and west brightness profiles for selected 1985 emission data and image cloud data. The spatial profiles both east and west of lo 
for the sodium D 2 emission brightness in units of kilorayleighs are shown as a function of the distance along the observing slit from the center of 
lo. Five emission observations identified by their satellite geocentric phase angles are shown by the different symbols. These five profiles occur 
when lo is east of Jupiter and past the satellite phase angle where the forward cloud has its symmetric turning point so that the trailing cloud 
profiles are all to the east (left) of lo and the forward cloud profiles are all to the west (right) of lo. A power law fit to each profile is also shown. 
At larger distances from lo, an envelope for the east-west D 2 emission profile acquired from sodium image data is shown by the shaded area. For 
the trailing profile, the shaded area is divided into two parts: the lower area corresponding to sodium cloud data when the directional feature is 
oriented either north or south and the upper area corresponding to the directional feature oriented along the east-west direction (i.e., the 
null condition). 


priate bounds for the brightness profiles when lo is some- 
what near the elongation points of its orbit. For the trailing 
cloud, two different shaded areas are shown in Fig. 2 for 
the two different basic orientations of the directional fea- 
ture: (1) lower area, when the directional feature is inclined 
either north or south, and (2) upper area, when the direc- 
tional feature is at or near the null locations. As expected, 
the shaded area for the directional feature near the null 
locations is both brighter and less steep than the shaded 
area for the directional feature with either a significant 
north inclination or a significant south inclination. At larger 
distances from lo (>30 i?i Q ), however, note that both 
shaded areas for the trailing sodium cloud are dimmer and 
more closely confined to lo than the shaded area for the 
forward cloud. In the forward cloud, all the emission 
brightness profiles are fairly tightly confined and have a 
slightly steeper slope than the shaded area due to increas- 
ing lo geometric plase angle. In the trailing cloud, the 
brightest and least steep of these profiles is for the emission 
4 (lo phase angle 87.7°) acquired for the directional feature 
at the null condition, while the next brightest profile is 
for the emission 5 (lo phase angle 100.6°), acquired —1.5 
hr later. 

4. ANALYSIS OF THE OBSERVATIONS 


taken. Collectively, the eclipse measurements for the co- 
rona near lo, the emission measurements that extend into 
the near sodium cloud, and the sodium cloud image derived 
profiles that reach to distances of ±100 R io , provide a set 
of spatially overlapping observations that will be used to 
study and constrain the initial velocity dispersion of the 
sodium source atoms at the exobase. In the modeling analy- 
sis, one-dimensional profiles are calculated using the nu- 
merical sodium cloud model of Smyth and Combi 
(1988a, b), where the electron impact ionization sink for 
sodium is determined for a 7° tilted corotating plasma 
torus with an offset-dipole planetary magnetic field in the 
presence of a nominal (i.e., —2.8 mV m _1 in Io’s frame) 
east-west electric field. A System III longitudinal asymme- 
try, although present in the torus ion emission, is not in- 
cluded but deferred to a later time when the electron de- 
pendence is available. 

To investigate the nature of the initial velocity dispersion 
of the sodium source, two different source flux speed distri- 
butions discussed earlier by Smyth and Combi (1988b; see 
their Appendix D) are considered: (1) a Maxwell- 
13 oltzmann flux distribution and (2) a modified-sputtering 
flux distribution. The Maxwell-Boltzmann flux distribu- 
tion <l>(v ; T) is based on the Maxwell-Boltzmann velocity 
distribution and is defined as 


Modeling analysis of the one-dimensional sodium distri- 
bution described in the previous section will now be under- 


<f)(v, T ) = 4 > o 


*y — f— y e -^ iv T> 2 

ReJ Vt \VtJ 


( 1 ) 
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Speed (km/s) 

FIG. 3. Flux speed distribution functions for sodium at Io’s exobase. Maxweli-Boitzmann flux speed distributions for sodium are shown for 
most probable speeds, u m , of 1.3 km sec -1 and 2.0 km sec -1 . Modified sputtering flux speed distributions are also shown for a = 3 and a most 
probable speed of 1.0 km sec -1 , for a = 7/3 and a most probable speed of 0.5 km sec -1 , and for a = 2 and a most probable speed of 0.4 km sec' 1 . 
All of the flux speed distributions are normalized to unit area under the curve. 


where u T = V2kT/m is the most probable speed of the 
velocity distribution for an atom of mass m. The Maxwell- 
Boltzmann flux distribution is proportional to the local 
velocity integrated flux referenced here to the satellite 
radius R s not the exobase radius R B and depends upon 
one parameter, the exobase temperature T (or alterna- 
tively U j), whic h determines both the most probable speed 
v m — V3 kT/m and the speed dispersion of the flux distribu- 
tion. The modified-sputtering flux distribution <f>(v; a, y b , 
u M ) is proportional to the local velocity integrated flux 
and depends upon three parameters, an exponent a and 
two velocity parameters v b and u M > 




( 2 ) 


where D(a , i/ M /u b ) is a normalization constant (see Smyth 
and Combi 1988b). The exponent a primarily determines 
the dispersion of the distribution, which has a greater high- 
speed population as a decreases. The exponent a has a 
value of 3 for a classical sputtering distribution (i.e., a 


complete collisional cascade process) and a value of 7/3 
for a Thomas-Fermi modified-sputtering flux distribution 
(i.e., the limit of a single elastic collisional ejection process), 
where the latter distribution is based upon a Thomas- 
Fermi differential scattering cross section. The velocity 
parameter v b is related nonlinearly to the most probable 
speed v m of the flux speed distribution and primarily deter- 
mines v m (see Smyth and Combi, 1988b, Appendix D). The 
velocity parameter primarily determines the maximum 
speed for the flux distribution and depends upon the maxi- 
mum relative speed (and masses) of the plasma torus ion 
and sodium atom. For different values of their parameters, 
two Maxwell-Boltzmann flux distributions and three mod- 
ified-sputtering flux distributions are shown in Fig. 3 and 
will be utilized in the subsequent modeling analysis. 

In calculating the column density and the D 2 emission 
brightness in the numerical sodium cloud model, a smaller 
two-dimensional sky-plane grid centered on Io (±15 Ri 0 ) 
is used to cover a spatial scale near the satellite more 
appropriate to the eclipse data while a much larger two- 
dimensional sky-plane grid centered on Io is used to cover 
a larger spatial scale more appropriate for the emission 
data and the sodium cloud image data. A one-dimensional 
profile for the eclipse data is obtained from the smaller 
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FIG. 4. Model calculations for the Io eclipse data using a Maxwell-Boltzmann flux speed distribution. The atomic sodium column density 
profile near Io determined from the 1985 eclipse data by Schneider et al (1991) is shown by the open circles. The modei calculated column density 
profiles are shown by solid dots for the (cylindrically averaged) corona, by solid triangles for the forward cloud along the east-west slit direction, 
and by solid squares for the trailing cloud along the east-west slit direction. These column density profiles were calculated using the Io sodium 
cloud model of Smyth and Combi (1988b) for their case C description of the plasma torus and for an Io geocentric phase angle of 92.9° and an Io 
System III longitude angle of 48.6°, which are similar to the emission 4 observation conditions in Table II. Sodium was ejected uniformly from an 
assumed exobase of 2600 km radius with a velocity dispersion for a Maxwell-Boltzmann flux distribution, where in (a) u m = 1.3 km sec -1 and 
<f> 0 = 3.0 X 10 8 atom cm" 2 sec -1 , and in (b) v m = 2.0 km sec -1 and = 1.8 X 10 8 atom cm -2 sec -1 (see text). 


two-dimensional sky-plane grid by extracting an average 
radial profile. This average radial profile (called the calcu- 
lated eclipse profile) will be denoted by the filled circles 
in Figs. 4-7. A one-dimensional east-west D 2 brightness 
profile (and also a corresponding column density profile) 
for the emission data and the sodium cloud image data is 
obtained from the larger two-dimensional sky-plane grid 
by selecting only the east-west grid elements that occur 
in the grid row containing Io. In Figs. 4-7, the calculated 
east -west brightness and column density profiles are de- 
noted by filled triangles for the forward cloud profile and 
by filled squares for the trailing cloud profile. To construct 
an eclipse or east-west profile, monoenergetic model cal- 
culations are performed for 18 different nonuniformly 
spaced speeds ranging from 0.4 to 10 km sec" 1 . Profiles for 
speeds beyond 10 km sec" 1 are determined by an inverse 
speed extrapolation of the model results. The individual 
profiles for the different speeds are appropriately weighted 
for a given source flux speed distribution and then added 
to obtain the final profile. Model calculations are per- 
formed for an Io geocentric phase angle of 92.9° and an 
Io System III longitude angle of 48.6°. These satellite condi- 
tions are similar to those for the emission 4 and eclipse 2 


observations of Table II, which are the observations closest 
to the eastern elongation point. This choice is also appro- 
priate for all the eclipse data, which has no discernible 
dependence on these two Io related angles, and for the Io 
sodium cloud image data which have east-west profile 
areas in Fig. 2 that are representative of the satellite near 
its orbital elongation points. Modeling analysis results are 
summarized in Table V and discussed below. 

For the first Maxwell-Boltzmann flux distribution in Fig. 
3 with a most probable speed of v m = 1.3 km sec" 1 (i.e., 
an exobase temperature of —1560 K) and with a flux <£ 0 
of 3.0 X 10 8 atoms cm" 2 sec" 1 (i.e., a total source of 
—1.2 X 10 26 atoms sec" 1 ), the model calculated eclipse 
profile (filled circles) in Fig. 4a provides an excellent fit 
within the Lagrange sphere to the eclipse observations 
(open circles) and also compares very favorably with the 
east-west column density profiles calculated for the for- 
ward (filled triangle) and trailing cloud (filled squares). 
This fit verifies and is similar to the earlier 1500 K Maxwell- 
Boltzmann flux distribution fit of Schneider et al (1991) 
noted in Section 2. Beyond the Lagrange sphere in Fig. 
4a, however, all three of these calculated profiles fall below 
the eclipse observations, which is considered less accurate 
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RADIAL DISTANCE FROM 10 (R, 0 ) 



RADIAL DISTANCE FROM lO (R j0 ) 

FIG. 5. Model calculations for the east-west D 2 brightness profiles using a Maxwell-Boltzmann flux speed distribution. The east-west D 2 
brightness profiles near Io in both the trailing and forward cloud directions as determined by the emission 4 data of Schneider et ai (1991) are 
shown by the open circles. The east-west profile envelopes in both the trailing and forward cloud directions as determined from the sodium cloud 
image data are shown by the shaded areas (see Fig. 2 caption). The descriptions for the calculated profile symbols, the sodium cloud model and 
plasma torus, and the Maxwell-Boltzmann flux distribution in (a) and (b) are the same as in the caption of Fig. 4. 
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RADIAL DISTANCE FROM IO (R j0 ) 


FIG. 6. Model calculations for the eclipse data using a modified sputtering flux speed distribution. The atomic sodium column density profile 
near Io determined from the 1985 eclipse data by Schneider et ai (1991) is shown by the open circles. The model calculated column density profiles 
are shown by solid dots for the (cylindrically averaged) corona, by solid triangles for the forward cloud along the east-west direction, and by solid 
squares for the trailing cloud along the east-west direction. These column density profiles were calculated using the Io sodium cloud model of 
Smyth and Combi (1988b) for their case C description of the plasma torus and for an Io geocentric phase angle of 92.9° and an Io System III 
longitude angle of 48.6°, which are similar to the emission 4 observation conditions in Table II. Sodium was ejected uniformly from an assumed 
exobase of 2600 km radius with a velocity dispersion for a modified sputtering flux distribution, where in (a) a = 3, v m = 1.0 km sec' 1 , and <&,= 
3.2 x 10 8 atom cm -2 sec -1 , in (b) a = 7/3, v m = 0.5 km sec -1 , and <&, — 4.2 X 10 8 atom cm -2 sec' 1 , and in (c) a — 2, v m = 0.4 km sec' 1 , and = 
4.7 X 10 8 atom cm' 2 sec' 1 . 
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at these distances. At and beyond about 8 R lo , the calcu- 
lated east-west forward (filled triangle) and trailing (filled 
squares) profiles rise above the calculated eclipse profile 
(filled circles) because the column density is no longer 
spherically symmetric about Io, with the forward cloud 
profile having the largest column density and showing a 
distinct change in its slope compared to the trailing cloud 
profile. The corresponding model profiles for the D 2 emis- 
sion brightness are given in Fig. 5a. For both the forward 
and trailing profiles, the calculated eclipse and calculated 
east-west profiles are in good agreement with each other 
inside the Lagrange radius, with a maximum brightness of 
about 200 kR near the exobase. The calculated east -west 
profile threads the three emission 4 data points for the 
forward cloud, but falls well below the emission 4 data 
points in the trailing cloud. For both the forward and trail- 
ing clouds at larger radial distances, the calculated east- 
west profiles fall well below the areas for both the forward 
and trailing cloud images. This behavior indicates that 
there is a large deficiency in the high-speed population for 
this source flux speed distribution. 

Model calculations were therefore performed for the 


second Maxweil-Boltzmann flux distribution in Fig. 3 with 
a higher most probable speed of v m = 2.0 km sec -1 (i.e., 
an exobase temperature of —3690 K) and with a flux 
4>o of 1.8 x 10 s atoms cm -2 sec" 1 (i.e., a total source of 
—0.75 X 10 26 atoms sec" 1 ) and are shown in Figure 4b and 
Figure 5b. For the D 2 emission brightness profiles in Figure 
5b, the calculated east- west profile now threads the center 
of the forward cloud image area for a radial distance up 
to about 70 Rio and the lower trailing cloud image area 
for a radial distance of about 25 R l0 before it falls off 
too steeply. This improved fit at larger radial distances, 
however, reduces the D 2 emission brightness at the exobase 
to about 80 kR in Figure 5b and causes the calculated 
eclipse profile in Figure 4b to fall below the measured 
eclipse profile for radial distances inside about 3 R Io . The 
Maxwell-Boltzmann flux distribution therefore cannot fit 
both the corona profile near Io and the sodium cloud east- 
west profiles at large distances from the satellite. A flux 
distribution that has a broader dispersion with enhanced 
populations for both the low-speed and high-speed atoms 
is required. The three modified-sputtering flux distribu- 
tions in Fig. 3, which have a broader dispersion, are thus 


TRAILING CLOUD FORWARD CLOUD 




RADIAL DISTANCE FROM IO (R !0 ) 


FIG. 7. Model calculations for the east-west D 2 brightness profiles using a modified sputtering flux speed distribution. The east-west D 2 
brightness profiles near Io in both the trailing and forward cloud directions determined from the emission 4 data of Schneider et ai (1991) are 
shown by the open circles. The east-west profile envelopes determined from the sodium cloud image data are shown by the shaded areas (see 
caption of Fig. 2). The descriptions for the calculated profile symbols, the sodium cloud model and plasma torus, and the modified sputtering flux 
distribution in (a), (b), and (c) are the same as in the caption of Fig. 6. 
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TABLE V 

Summary Comparison of Modeled and Observed East- West Sodium Profiles for Different Flux Speed 

Distributions at lo’s Exobase 



Corona 

Forward Cloud 

Trailing Cloud (not null) 

Trailing Cloud (null) 

Observations: E/W Radial Interval (Rio): 

1.4-6 

6-10 

10- too 

6-10 

10- too 

6-10 

10-100 

Distribution Exobase 

Speed Peak Source Rate 

(km/s) (10« atoms/s) 

Fits 

Eclipse 

Column 

Profile 

Fits 

Near lo 

Emission 

Profile 

Fits 
Far 
Cloud 
E/W Profile 

Fits 

Near Io 

Emission 

Profile 

Fits 
Far 
Cloud 
E/W Profile 

Fits 

Near Io 

Emission 

Profile 

Fits 
Far 
Cloud 
E/W Profile 

1. Maxwell Boltzmann Flux Distribution 








1.3 1.24 

YES 

YES 

too low 

little low 

too low 

too Low 

too low 

2.0 0.75 

too low 

YES 

slightly low 

YES 

little low 

little low 

too low 


2. Coliisional Cascade Flux Distribution 
(classical sputtering) 

a = 3 1.0 1.32 YES YES tiny low YES tiny low little low too low 


(incomplete cascade: higher velocity tail) 


a = 7/3 0.5 

1.74 

YES 

YES 

YES 

YES 

YES 

little low 

too low 

a - 2 0.4 

1.90 

tiny low 

too high 

too high 

too high 

too high 

YES 

YES 


considered in the remainder of the paper with model calcu- 
lations presented in Figs. 6 and 7. 

Model calculations for a classical sputtering flux distribu- 
tion (a = 3) and a modified-sputtering flux distribution 
(a = 7/3) are presented in Figs. 6a and 6b for the eclipse 
observations and in Figs. 7a and 7b for the east-west D 2 
emission brightness profiles. For these two flux distribu- 
tions, the most probable speeds are, respectively, 1.0 km 
sec” 1 and 0.5 km sec” 1 , and the sodium fluxes <£ 0 are, respec- 
tively, 3.2 X 10 8 atoms cm" 2 sec" 1 (i.e., a total source of 
-1.3 X 10 26 atoms sec” 1 ) and 4.2 X 10 8 atoms cm' 2 sec" 1 
(i.e., a total source of —1.7 X 10 26 atoms sec" 1 ). From the 
exobase to radial distances of —8 R\ 0 , just beyond the 
Lagrange radius, both sputtering flux distributions provide 
a very good fit in Figs. 6a and 6b to the observed eclipse 
column density profile (open circles) and correspond to 
an exobase D 2 emission brightness of about 150 kR in Figs. 
7a and 7b. For the classical sputtering flux distribution in 
Fig. 7a, the calculated D 2 emission brightness profile for 
the forward profile is slightly above the measured data 
point (open circles) inside the Lagrange radius, matches 
the two measured data points beyond the Lagrange radius, 
and then threads the forward cloud image area nicely be- 
tween about 20 R lo and 80 R lo before it falls too rapidly 
and drops below this area. An excellent fit for the forward 
profile is, however, provided by the modified sputtering 
distribution ( a - 7/3) in Fig. 7b where the calculated D 2 
emission brightness profile matches the measured data 


points (open circles) both inside and beyond the Lagrange 
radius as well as nicely threading the forward cloud image 
area all the way to 100 # Io . For the trailing cloud, the 
calculated D 2 emission brightness profile for the classical 
sputtering flux distribution in Fig. 7a matches the measured 
data point inside the Lagrange radius, is slightly below the 
two measured data points outside the Lagrange radius, 
and then threads the lower of the two trailing cloud image 
areas nicely between about 15 Ri Q and 35 R io before it falls 
too rapidly and drops below this area. An excellent fit for 
the trailing profile is, however, provided by the modified 
sputtering distribution (a = 7/3) in Fig. 7b where the calcu- 
lated D 2 emission brightness matches the measured data 
point inside the Lagrange radius, is slightly below the two 
measured data points outside the Lagrange radius, and 
then threads the lower (non-null) trailing cloud image area 
nicely all the way to 100 Ri 0 . It is particularly noteworthy 
that the isotropic ejection of sodium from the exobase 
with a modified sputtering flux distribution with a = 7/3 
provides a complete fit to the combined eclipse, emission, 
and forward/trailing sodium cloud image profile data for 
this non-null condition from 1.4 to 100 /? Io . 

In order to fit the trailing cloud (upper area) profile for 
the directional feature at the null condition, it is then clear 
that a flux distribution is required with an even more en- 
hanced higher-speed population (—20 km sec" 1 ) than the 
modified sputtering flux distribution with a = 7/3. Since 
the modified sputtering flux distribution for a = 7/3 corre- 
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sponds to the limit of a single collision cascade process 
described by a Thomas-Fermi cross section (see Smyth 
and Combi 1988b), reducing the value of a to a smaller 
value becomes somewhat physically questionable but will 
be used here for the purposes of simply illustrating the 
impact of a more enhanced higher-speed sodium popula- 
tion in the model calculation. As discussed earlier, this 
higher-speed sodium is thought to be nonisotropically 
ejected from Io’s exobase and attributed to some combina- 
tion of direct collisional and lower-velocity charge ex- 
change ejection. Choosing the modified sputtering flux dis- 
tribution with a — 2 in Fig. 3 which has a most probable 
speed of 0.4 km sec -1 and selecting an isotropic exobase 
source rate of 1.9 X 10 26 atoms sec -1 (i.e., a flux <f>o of 
4.7 X 10 8 atoms cm -2 sec -1 ), the model-data comparison 
is shown in Fig. 6c for the eclipse column density and 
in Fig. 7c for the east-west D 2 emission brightness. The 
sputtering flux distribution provides a reasonably good fit 
to the observed column density data points in Fig. 6c with 
only a small departure very near the exobase and produces 
a column density profile beyond 10 R io that is significantly 
enhanced compared to the a = 7/3 case in Fig. 6b. In-Fig. 
7c, this enhancement in the forward cloud is obvious, where 
the calculated D 2 emission brightness profile is significantly 
above the measured data points both inside and outside 
the Lagrange radius and is above or in the very top of the 
forward cloud image area all the way to 100 R io . The 
additional enhanced high-speed population of the a = 
2 modified sputtering flux distribution is too large and 
therefore not consistent with the observed forward profile. 
In contrast for the trailing cloud in Fig. 7c, the calculated 
D 2 emission brightness profile matches the measured data 
points inside and outside of the Lagrange radius very well 
and then threads the upper of the two trailing cloud image 
areas nicely all the way to ~90 R lo . This demonstrates that 
the trailing cloud can be fitted with an enhanced higher- 
speed population of sodium atoms in the flux distribution. 
It also immediately demonstrates that the flux distribution 
at the exobase must be nonisotropic with the enhanced 
high-speed population weighted toward vector directions 
that will preferentially populate the trailing cloud. As dis- 
cussed in Section 2, this nonisotropic requirement for a 
flux distribution for speeds of ~20 km sec" 1 is similar to 
the conclusions reached by earlier modeling analyses (Pil- 
cher etal 1984; Smyth and Combi 1991; Wilson and Schnei- 
der 1995). 

5. DISCUSSION AND CONCLUSIONS 

The composite spatial information for sodium obtained 
by combining the eclipse observations (radial distances 
from Io of 1.4 to ~10 i? Io ), the emission observations (east- 
west distances of ±4 to ±30-40 7? Io ), and the sodium cloud 
observations (east- west distances of ±10 to ±100 Z? Io ) has 



FIG. 8. Two-dimensional nature of the sodium column density near 
Io. Contours for the two-dimensional column density near Io are shown 
in the sky-plane of the earth as determined from the sodium cloud model 
calculation for the modified sputtering flux speed distribution described 
in Fig. 6(b) for a - 7/3. The vertical and horizontal directions are the 
projected directions that are, respectively, perpendicular and parallel to 
the semimajor axis of the Io’s orbital ellipse on the sky plane. Io’s location 
and size are shown to scale by the black circle. The sodium column 
density contours in units of 10 11 atoms cm -2 are, from inside to outside, 
5, 2, 1, 0.5, 0.2, 0.1, and 0.05. 


been analyzed to extract a basic description for the flux 
speed distribution at the satellite’s exobase. An isotropic 
modified-sputtering flux speed distribution in Fig. 3 with 
a = 7/3, a most probable speed of 0.5 km sec -1 , and a 
source strength of 1.7 X 10 26 atoms sec" 1 provides a very 
good fit to these composite observations when the direc- 
tional feature is either north or south and hence not con- 
tributing to the east-west profile of the trailing cloud. It 
is remarkable that these observations, acquired by a num- 
ber of ground-based programs over very different spatial 
scales and at different times during the 1976-1985 decade, 
are so self-consistent. Near Io, the two-dimensional sodium 
column density produced by this modified sputtering distri- 
bution as calculated by the sodium cloud model in the 
profile analysis above is shown in Fig. 8 and can be seen 
at larger distances from Io to become nonspherical and 
more confined near the satellite plane. This flattening near 
the satellite plane is the merging of the near Io corona 
into the sodium cloud and is caused naturally by orbital 
dynamics beyond the satellite Lagrange sphere where the 
gravity of Jupiter is dominant. The forward cloud portion 
of the east-west emission data profiles has a rather tightly 
confined slope that, in the absence of the trailing cloud 
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FIG. 9. Total source rate speed distribution function for sodium at Io’s exobase. The total source rate speed distribution function at Io’s exobase, 
in units of 10 26 atoms sec -1 (km/sec)" 1 , is composed of three separate source rate speed distributions as discussed in the text and is shown for two 
different source strengths for the higher-speed zenocorona source centered about 57 km sec -1 . The lower (solid line) and upper (dashed-dot line) 
curves correspond, respectively, to the sodium zenocorona higher-speed source rates of 2.2 x 10 26 atoms sec -1 (Smyth and Combi 1991) and 4 X 
10 26 atoms sec' 1 . The two source rates for the higher-speed zenocorona source are shown to exhibit its typical time-variable source strength range 
of ~2-4 X 10 26 atoms sec' 1 as reported by Flynn et al (1994). The decomposition of the solid curve into its three separate source rate speed 
distributions is shown in the cutout and is determined by combining (1) the isotropic modified sputtering source rate distribution (dotted line in 
the cutout) for a = 7/3, v m = 0.5 km sec' 1 and a source strength of 1.7 X 10 26 atom sec' 1 , (2) the nonisotropic lower-speed source rate distribution 
(short dashed line in the cutout) for the sodium zenocorona and directional feature centered about 20 km sec' 1 , with a source strength of 1.1 X 
10 26 atoms sec' 1 as determined by Smyth and Combi (1991), and (3) the nonisotropic higher-speed source rate distribution (longer dashed line in 
the cutout) for the sodium zenocorona centered about 57 km sec -1 , with a charge exchange source strength of 2.2 X 10 26 atoms sec' 1 as determined 
by Smyth and Combi (1991). 


enhancement at the null condition, is less steep and is 
brighter than the trailing cloud profiles. In order, however, 
to reproduce the extended east-west profile in the trailing 
sodium cloud when the directional feature is in the satellite 
plane (i.e., the null location), additional nonisotropic high- 
speed sodium is required. 

The sodium atoms ejected from Io’s exobase as de- 
scribed above by the modified sputtering flux distribution 
have speeds primarily in the range from 0 to a few tens 
of km sec -1 . This neutral flux distribution represents the 
spatially integrated effect of the incomplete collisional cas- 
cade process that occurs from the collisional interactions 


of heavy ions in the corotating plasma torus with neutrals 
in Io’s atmosphere. This flux speed distribution can be 
alternatively described as a source rate speed distribution 
by multiplying it by the satellite surface area. In addition to 
these ion-neutral elastic collisional encounters, resonance 
charge exchange between plasma torus sodium ions and 
neutral sodium in Io’s atmosphere is also responsible for 
producing a sodium source with higher speeds relative to 
Io. These speeds are centered about the corotational ion 
speed (—57 km sec” 1 ) relative to Io’s motion and have a 
dispersion reaching from several tens of km sec” 1 to —100 
km sec” 1 . Such high speed sodium (<80 km sec” 1 ) has 
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recently been observed near Io by Cremonese et al (1992). 
As discussed in Section 2, this higher-speed nonisotropic 
source of sodium together with the lower speed (—15-20 
km sec -1 ) nonisotropic source for the directional feature 
form the source for the sodium zenocorona or magneto- 
nebula. Earlier modeling studies (Smyth and Combi 1991; 
Flynn et al 1992) indicated that the higher-speed source 
was —2 X 10 26 atoms sec -1 while the lower speed source 
was —1 X 10 26 atoms sec -1 . More recent observations and 
analysis (Flynn et al 1994) have shown that the higher- 
speed sodium source is time variable and in the range 
—2-4 X 10 26 atoms sec -1 . A total source rate speed distri- 
bution for sodium at Io’s exobase has hence been con- 
structed by combining the modified sputtering source rate 
distribution determined in this paper with the two source 
rate distribution for the zenocorona as given by Smyth and 
Combi (1991) and is shown in Fig. 9. The lower (solid line) 
and upper (dashed-dot line) curves correspond, respec- 
tively, to the sodium zenocorona higher-speed source rates 
of 2.2 X 10 26 and 4 X 10 26 atoms sec -1 . Total source rate 
speed distribution functions at Io’s exobase expected for 
other atomic species, such as K, O, and S, can be -con- 
structed in a similar fashion by adopting the estimated 
source rates given by Smyth and Combi (1991). 

Future studies for the sodium flux speed distribution at 
Io’s exobase are anticipated using much larger data sets 
now available for east-west and north-south sodium emis- 
sion observations. It will then be possible to analyze the 
combined spatial and spectral information and refine the 
nonisotropic nature of the flux distribution and also to 
search for other possible east-west and System III related 
variations. Once this information is determined for sodium, 
the implications for the more abundant species in Io’s 
atmosphere will be particularly important in other related 
studies for the many-faceted and complex phenomena in 
the Io-Jupiter system. 
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